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The existence of gauge (Gribov) copies disturbs the usual Faddeev-Popov quantization procedure 
in the Landau gauge. It is a very hard job to treat these in the continuum, even in a partial manner. 
A decent way to do so was worked out by Gribov, and later on by Zwanziger. The final point 
was a renormalizable action (the Gribov-Zwanziger action), implementing the restriction of the 
path integration to the so-called Gribov region, which is free of a subset of gauge copies, but 
not of all copies. Till recently, everybody agreed upon the fact that the restriction to the Gribov 
region implied a infrared enhanced ghost, and vanishing zero momentum gluon propagator. We 
discuss how the Gribov-Zwanziger action naturally leads to the existence of vacuum condensates 
of dimension two. As it is very common, such condensates can seriously alter the dynamics. In 
particular, the Gribov-Zwanziger condensates give rise to a gluon propagator with a finite but non- 
vanishing zero momentum limit, and reconstitute a nonenhanced ghost. We call this the refined 
Gribov-Zwanziger framework. The predictions are in qualitative agreement with most recent 
lattice simulations, and certain solutions of the Schwinger-Dyson equations. A crucial feature 
of the Gribov-Zwanziger framework is the soft (controllable) breaking of the BRST symmetry. 
We also point out that imposing the Kugo-Ojima confinement criterion on the Faddeev-Popov 
theory as a boundary condition from the beginning leads to the same partition function as of 
Gribov-Zwanziger, with associated BRST symmetry breaking. This clouds the interpretation of 
the Kugo-Ojima criterion in se. 
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1. The Gribov-Zwanziger framework 

As it is well-known, the classical SU (N) Yang-Mills action, 

Sym = ^Jd 4 xF« v F« v , (1.1) 
displays an enormous invariance w.r.t. 

A^A^+Dfco 17 (1.2) 

with co a arbitrary (infinitesimal) functions 1 . In order to quantize this gauge theory, one needs in 
principle to select a single gauge representant A^ by fixing the gauge. A popular choice is to restrict 
to the transverse degrees of freedom by imposing the Landau gauge, d^A"^ = 0. 

After some nice machinery at the level of the path integral, one ends up with the Faddeev- 
Popov gauge fixed action in the Landau gauge. 

Sym+fp = \J d 4 xF« v F£ v + f d 4 x (b a d^ +7*d^Dfc b ) , (1.3) 

which comes from the partition function Z by 

Z = J [dA] 8 (dA) detMe- SyM = J [dA] [db] [dc] [dc]e~ Sru+FP . (1 .4) 

Indeed, setting b a = iB a , one finds the functional analog of the Fourier representation of a 8- 
function. The accompanying Jacobian detM is lifted into the exponential using the ghost fields c a 
and c a . We introduced the Faddeev-Popov operator as 

M ah = -d^Df = (d^8 ab + g r%) . (1.5) 



Clearly, the local gauge invariance expressed by (|L2j) is lost. It is however replaced by the BRST 
invariance ^ |||], 

sAl = - {D,c) a , sc a = l -gf hc c h c c , sc" = b" , sb a = 0, 

with nilpotent generator s, s 2 = 0. This BRST symmetry is a cornerstone of (perturbative) gauge 
theories. 

In the above construction, it was silently assumed that there is only one solution on the 
gauge orbit of an arbitrary gauge field that obeys the Landau gauge, i.e. d^A^ = 0. Let us consider 
a gauge equivalent field A^, 

A^'=A^+Dfco b , (1.6) 
then this A^ could also obey the Landau gauge condition, d^A'^ = 0, if 

M ab (O b = 0, (1.7) 



We shall not speak about "large" gauge transformations which are not connected to the unity. 
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i.e. if the Faddeev-Popov operator has a zero mode. It was the insight of Gribov to show that such 
zero modes do exist, hence we encounter the problem of gauge copies [J3j] - This of course poses 
a serious problem for the Faddeev-Popov quantization procedure, where it was e.g. assumed that 
detM > 0, but since M has zero modes, it can switch sign 2 . For perturbation theory, these copies 
are pretty harmless, as discussed in e.g. [Q]. However, when one is interested in a region where 
nonperturbative physics could enter the game, one needs to take care of these copies. Therefore, it 
would be interesting to improve upon the original Faddeev-Popov procedure. 

This has been worked out by Gribov in his seminal paper in a steepest descent approximation. 
We quote here his result, which is an improved version of the Landau gauge fixed partition function, 
reading 

■ d <t,] e -SYM + FP+ffd 4 xAjzA_ (L8) 

This was achieved by writing the inverse Faddeev-Popov operator (with an external gauge field) as 

M-' = (-BD)-' = w -L m} (1.9) 

and inserting into the path integral a step function 6(1 — a(0,A)), which ensures that are are only 
integrating in A-space over configurations wherefore M is positive. The region 

Q = {A a }ll d ll A a il = Q,M ab >0} 

is known as the Gribov region, and it displays many nice properties, summarized in [H]. For 
example, Q. is convex, bounded in all directions in field space, and it has the crucial property that 
each gauge orbit intersects with Q.. 



After a few steps, one then arrives at (L8), whereby y 4 is a thermodynamic mass parameter, 
fixed by 

3 ,„ I A'q 1 



4 sN Jw?7Tf = L (U0) 

We notice here that this gap equation is ill-defined as the integral diverges. The condition 

ff(0,A) < 1 (1.11) 
is known as the Gribov no-pole condition. We notice here that if A becomes a quantum field 



and is integrated over, the expectation value of the operator ( |1.9[ ) corresponds exactly to the ghost 
propagator. 

By implementing the restriction to the Gribov horizon, one already avoids copies which are 
infinitesimally related to each other. Unfortunately, there are still other copies left inside Q., see 
e.g. A further restriction would be necessary, but it is unclear how this could be achieved in an 
analytical fashion. 

Having an improved partition function at our disposal, one has the ideal tool to study the 
theory. In particular, one can look at the elementary gluon and ghost propagator. For the gluon 



2 We recall that the 5-function is accompanied by the absolute value of its corresponding Jacobian. It is however 
unknown how to lift | detM into the exponential. 
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propagator, one finds at tree level 3 

<A«A«) p = S ab D(p 2 ) (V - EEj*^ , D(p 2 ) = -j^j, (1.12) 



while for the ghost, 

''' ' ' ' /7i oip 



(c a c b )=8 ab G(p 2 ), G(p 2 ) = ^ T i * (1.13) 



this at one loop. Implementing the gap equation leads to a(0) = 1, thus the ghost propagator 
exhibits an infrared enhancement. This should not come as a big surprise, since the gap equation 
( 1 . 1 Q[ ) is the quantum version of the no-pole condition ( |1.1 1| ). At p 2 = 0, we reach the boundary 
value a(0) = 1, which corresponds to the boundary (horizon) dQ. of the Gribov region £2. For the 
gluon propagator, we notice an infrared suppression, in particular D(0) = 0. 

So far, everything we have mentioned is restricted to the lowest (quadratic) order. After a tour- 
de-force, Zwanziger was able to extend Gribov's result to all orders [{7|, §J, giving the following 
action 

S h = Sym+fp + 7 4 J d 4 * h(x) , (1.14) 



with the horizon function 



h(x) = g 2 f abc A h l (M- l ) ad f dec A^. (1.15) 



The parameter y is now fixed by demanding that 

(h(x)) = d(N 2 -\), (1.16) 

the latter being known as the horizon condition. It can be easily checked that the Zwanziger re- 
sults reduce to those of Gribov at lowest order. We draw attention to the fact that the action is a 
nonlocal one. As nonlocal quantum field theories are hard to handle and interpret, and almost all 
available tools are intended for local quantum field theories, this looks troublesome. Fortunately, 
by introducing a suitable set of extra fields, we can reformulate ( 1.14 ) as 

Sgz = Sym+fp + Sh , (1.17) 

with 

5, = J d A x(l?p v (d v ^ c + gf bm A h v 9T') -<dv [d v co^ + gf hm A b v co; c ) 

-g (d v 7D%) f abm (D v c) b <p™ - fg (r bc A^(p bc +f abc A a ^ + ~ g (N 2 - 1) T 2 ) ) ,(1.18) 

which is known as the Gribov-Zwanziger action. The horizon condition ( |1.15| ) turns out to be 
equivalent to 

= 0<*(gf^Al(y + y)*)=2d{N 2 -\)f, (1.19) 



3 We set A 4 = 2g 2 Ny 4 
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where the artificial solution y 2 = needs to be discarded. We see that the horizon condition corre- 
sponds to a certain nontrivial d = 2 condensate. 

Having constructed a local action implementing the restriction to £2, one can really start to 
study the theory. In particular, we can mention that S GZ defines a renormalizable theory, which 
is a highly nontrivial issue [M 0, 10]. In particular, we draw attention to the fact that the factor 



in front of the pure vacuum term in ( 1.18| ) turns out to be uniquely fixed by the renormalization. 
If this would be not the case, as it is in general the case [O], it would be impossible to speak 
about the renormalized version of the horizon condition ( |1.19[ ), and thus of the gap equation ( l.lfjp . 
We have found this to be a truly remarkable feature of the Gribov-Zwanziger framework. For the 
propagators, one basically finds back the Gribov predictions at higher order too, i.e. an infrared 
suppressed gluon propagator and a enhanced ghost propagator. This can be proven using the Ward 



identities of the theory as done in [Bj], while it was also explicitly verified in works as [12J, 13Q. 



2. Fate of the BRST symmetry in the Gribov-Zwanziger framework 



Let us first mention that for 7=0, the physics described by the action (1.18) is completely 
equivalent to that of the original Faddeev-Popov theory. This can be appreciated since there exists 
a natural extension of the original BRST transformation to these fields, by setting 



CO] 



ab 



SCO] 



ab 











The action ( |1.18| ) is then an s-exact piece, Sh = s(. . .), while the extra fields come as ^-doublets, 



hence they are physically trivial Q14[]. In the case of interest, 7^0, and we arrive at a more peculiar 
conclusion. We find 



sS G z = g 



f J d 4 xf abc 



(D a ™c m ) 



9m +< 



(2.1) 



or we observe a soft BRST breaking. The breaking is called soft, as it is only quadratic in the 
fields. Therefore, it is perfectly under control at the quantum level, and one can even write down 
the corresponding softly broken Slavnov-Taylor identity [|D|], 



/ 



d 4 xi 



sr sr sr sr 

SK" SA" + ~8L' 



^ sr ia sr _ a sr a sr \ fA , 



(2.2) 



whereby [Ay • r] corresponds to insertion of the breaking operator Ay. A nice application is that 
this allows to prove that the gluon self energy is no longer transverse. This was explicitly checked 
to two loops in [|i"3|]. 

At the level of renormalizability, there is thus no problem associated with the softly broken 
BRST. However, as it is well know, the BRST also plays a crucial role in establishing the unitarity 



of gauge theories [16], including the definition of the physical states, which a fortiori must have a 
positive norm. Let us therefore first give a quick look at the usual Faddeev-Popov case. 



3. The Kugo-Ojima analysis 



In [16], building on work [17], a beautiful discussion of the unitarity issue based on the nilpo- 
tent BRST charge was given. When physical states \\jf p ) are defined as those states belonging to 
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the BRST cohomology 4 , \\\f p ) € co\\oyci{Qbrst) = ^Qbrst ' one can snow tnat on ty ^ e ^ trans- 
verse gluon polarizations are physical; the other two cancel with the (anti)ghost particles. This 
is a generic feature, mainly based on purely algebraic properties, stating that unphysical particles 
always decouple in quartets (the members of a quartet are related to each other by the BRST and 
ghost charge). Roughly speaking, this property remains preserved under time evolution, as the 
BRST charge commutes with the Hamiltonian as a symmetry generator. Said otherwise, it allows 
one to define a unitary ^-matrix. 

Kugo-Ojima continued their analysis, and also considered the global color charge, which can 
be written as 5 

Q a = J d\d i F^ + J d 3 x{ Qbrst i Dqc" } . (3.1) 

Assuming that gluons are infrared suppressed in the sense that the first volume-integral is well- 
defined, thus JfPxd'FZ = 0, and that / ^xIQersttDqc } is well defined too, we would have color 
confinement on physical states, due to Q BRST = 0. The big problem is of course to assure that the 
space integral / d 3 x{QsRST ,Doc a } is well-defined. It was discussed in [ |16| , 18] that this is the case 
if 

w(0) = -l, (3.2) 
where u(p 2 ) is the transverse part of a particular Green function, namely 

PuPv 



J dV> (D»c a {x)D v c h {<d)) = 8 ab U, v - P -^f\ u(p 



J " ' p 2 

An interesting relation now exists between the ghost propagator and this u(p ) in the specific case 
of the Landau gauge, first shown in [|TJ] and confirmed in [19], 

G ab (p 2 ) = (c a c b ) = (3.3) 
\ I p p 1 l + u(p z )+p z v(p l ) 

and recently rederived in a slightly different setting in the paper [p0|]. It is generally accepted that 
v(0) = [pot], and in a loop expansion, it is actually zero as one can learn from e.g. [|l2|], therefore 
one notices that u(0) = — 1 is realized when the ghost propagator is infrared enhanced, and vice 
versa. This explains the attention paid to the infrared behaviour of the ghost in the past, both from 
numerical and analytical viewpoint. We refer to e.g. [15] for an extensive list of literature. 

We can make a few comments. First of all, it was unclear at the time whether the Kugo-Ojima 
function u{k 2 ) could be renormalized. In the Landau gauge, such a proof was presented in [^l], p2| ]. 
However, it is clear that w(0) = — 1 is not realized in perturbation theory, hence one needs to resort 
to nonperturbative calculations tools. One such example are lattice gauge studies, and most recent 
results agree on the fact that m(0) = — 1 is not realized, either via direct studies or via the ghost 
propagator that appears to be not enhanced. Clearly, an important role in the analysis is played 
by the BRST charge, however it is unknown whether a well-defined BRST charge exists in the 
infrared. Moreover, the Kugo-Ojima analysis is purely based on the Faddeev-Popov action, and 
completely ignores the gauge copy problem. 

4 It is of no use to keep the BRST-exact states, Qbrst !■••)> as these are trivially annihilated by Qbrst> moreover, 
they have zero norm. 

5 We do not pay attention here to working in Minkowksi or Euclidean space. The ideas thereafter will be clear. 
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4. Combining Gribov-Zwanziger with Kugo-Ojima 

Having summarized in a nutshell the Gribov-Zwanziger approach and Kugo-Ojima results, it 
is instructive to compare both issues a bit. As already mentioned, the Gribov restriction suppresses 



the gluon in the infrared, which should make the first part in (3. 1 ) well-defined. Simultaneously, the 
ghost is enhanced when the horizon condition is implemented, so apparently the Gribov-Zwanziger 
framework fulfills all the Kugo-Ojima requirements to ensure confinement. The correspondence 



becomes even more apparent when we compare (1.13) and (3.3), giving that a(0) = — u(0). The 



no-pole condition translates into a(0) = 1 at the quantum level, hence u(0) = —1. Recently, the 
correspondence between the nonlocal horizon function ( 1.14 ) and the "pole function" o(k 2 ) was 
also shown at higher orders p3j\. 

The link between u(p 2 ) and the horizon condition can also be formalized. Using the Ward 
identities of the Gribov-Zwanziger action, one shows that [||] 

gf^Alto + lfl*) =2fJ d 4 x(D^(x)D^(0)) (4.1) 
= -2f(N 2 -\)((d-l)u(0)-\). (4.2) 



Hence, implementing the horizon condition (1.19) is equivalent with setting u(0) = — 1. 



We notice that the horizon condition can thus be connected to the correlator 

(D fl c a (x)D^c a (y)), (4.3) 

which can be done in a renormalizable setting. Other versions of the horizon function can be found 
in the literature but to our understanding this would correspond to a connection with the 
correlator 

\gf ahc A'^\x)){ g f^A p ^){y)) (4.4) 

which is however not renormalizable [24]. We claim that the requirement of renormalizability fixes 
the "choice" in the horizon function in a unique way. 

From the above considerations, it appears that there exists something like a GZKO framework, 
where everything fits wonderfully together. However, there is a big hole in this reasoning. As 
we already explained, the Gribov-Zwanziger restriction breaks the BRST, while the Kugo-Ojima 
analysis indispensably needs the BRST. This means that actually, in the Gribov-Zwanziger theory, 
the meaning of the Kugo-Ojima criterion is unclear. 

A potential solution would be to construct a new BRST charge for the Gribov-Zwanziger 
theory. That this might be possible can be guessed from the fact that at the quadratic order, 

sS GZ = *(y 4 J d 4 ^^A^ oc /| d 4 xA,±d,c a , (4.5) 

and performing a partial integration, it appears that this breaking could be absorbed into a redefini- 
tion of the b a -field. A new BRST symmetry STor the complete theory was indeed found in 



but it turns out to be nonlocal, as already explained on general grounds in [15]. It remains un 



clear what a nonlocal BRST symmetry could be used for, it is even unclear what the corresponding 
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charge would be. Evidently, a new BRST would also need a reworking of the Kugo-Ojima analysis, 
if possible of course. 

Anyhow, till this moment, there is absolutely nothing proven about the nonperturbative unitar- 
ity of or confinement in the Yang-Mills theory in the Landau gauge, when relying on the Gribov- 
Zwanziger and/or Kugo-Ojima approach. All that one can say is that implementing the Gribov 
restriction leads to a violation of positivity in the spectral representation of the gluon propagator 
[15]. This violation traces back the presence of complex poles in the propagator ( 1.12 ), and as such 
the gluon cannot describe a physical particle. This violation of positivity has also been found in 
other approaches, like in L28h. 



5. Imposing Kugo-Ojima as a boundary condition 

In this section, we wish to implement the Kugo-Ojima criterion as a boundary condition in 
the theory. This was worked out in p9|], and it was inspired by works in which this is also done 
to select a specific solution. A very powerful tool to extract nonperturbative information from a 
quantum field theory is by solving their quantum equations of motion, being the Schwinger-Dyson 
equations. Although this is a horrendous task, much progress has been made in recent years, we 
refer to e. g. [fHj, || ||] for an overview of the literature and results. There exists whole classes 



of solutions for these equations, and one must impose certain conditions to discriminate between 
them. A possibility is to impose a condition on the ghost propagator in the deep infrared. One can 



imagine imposing an infrared enhanced ghost [31 1, which corresponds to the a = case of 



[p 2 G(p 2 )] p i = a>0. (5.1) 
It has been speculated that this choice of a can be related to different nonperturbative implementa- 



tions of the Landau gauge [[33J . 

Since imposing a constraint on a propagator in an interacting quantum field theory corresponds 
to demanding something highly nontrivial at the quantum level, one can imagine that such a con- 
straint might corrupt e.g. some symmetry properties of the original action. We have explained in 



[|29p that for a consistent treatment, one would need to impose the constraint from the start, in 
order to have a decent way to discuss the symmetries of the constrained theory. It turns out that 
this is possible in case of the Kugo-Ojima criterion. More precisely, it is possible to impose that 
m(0) = — 1, or thus that the ghost is infrared enhanced, and this in a way that is stable under quan- 
tum corrections. The latter is an important ingredient, as if this would not be possible, it would 
make no sense trying to enforce the theory's behaviour in such a way. In practice, we stalled from 
the following rather unconventional form of the Faddeev-Popov action, 

S'ym + fp = S ym+ fp + J d 4 *(V-d v (d v cp£ + gf bm A b v <p™) -a%d v (d v co"; + gf hm A h v co™) 

-g(d v W a <)f abm (D v c) b <p™), (5.2) 



which is completely equivalent with the original one, (1.3) (see Section 1). Next, we introduce a 
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multiplier y 2 as follows, 

S'ym+fp ~^ Sko = S'ym+fp 

- |d 4 x( 7 2 ^(r^(<+^) + l(^ 2 -l)7 2 )), (5.3) 

and impose the condition = 0, which enforces that 

' f*c A a (<p + <p)J<A = 2 ^ 2 - 1)/. (5.4) 



As we have already discussed before, this is in return sufficient to ensure the ghost enhancement, 
through u(0) = — 1. 

The reader shall of course observe the similarity between this result and the one concerning 
the Gribov-Zwanziger action. We want to stress here that currently, we are not saying anything 
yet about copies or treatment thereof 6 . What we are doing is imposing the constraint ( |5~4| ) on top 
of the normal Faddeev-Popov action. We then obtain an action that automatically implements the 
Kugo-Ojima criterion as a boundary condition, in way that is consistent at the quantum level, due 



to the renormalizability of the associated action (5.3). We afterwards notice that the eventual action 



is indeed equivalent to the one constructed by Gribov and Zwanziger. 

We conclude that imposing the Kugo-Ojima criterion on top of the Faddeev-Popov quantiza- 
tion softly breaks the BRST symmetry 7 . As such, the role of the condition w(0) = — 1 in explaining 
confinement becomes again difficult to understand. 

6. Exploring the Gribov-Zwanziger quantum dynamics 

In this concluding section, we shall spend a few words on the quantum dynamics of the Gribov- 
Zwanziger theory. The motivation arises from the comparison with lattice results, which essentially 
come to the agreement that the gluon propagator is suppressed in the infrared, but it does not vanish 



at zero momentum, while the ghost propagator does not seem to be enhanced [ 34 ] (see however 



e.g. [|35|, ^3, [37p for other views). 

It is perhaps important to notice that the restriction to the Gribov region £2 in fact corresponds 
to selecting the minima along the gauge orbit of the functional / d 4 xA 2 . The problem of residual 
copies is caused by the existence of multiple local minima. One should resort to finding the absolute 
minima In the continuum, it has been argued that it would be sufficient to restrict to Q. when 
calculating expectation values [|38|]. In the continuum, the restriction to CI is done via the Gribov- 



Zwanziger action ( |1.18| ). It is however not possible to simulate with the Gribov-Zwanziger action 
on a lattice. Numerically, one rather selects the gauge configuration that corresponds to the "best" 
minimum each time. It is always assumed that in the infinite volume limit, the lattice and the 
continuum version are essentially doing the same. 

This being said, in the past, using relatively small lattices, there seemed to be agreement 
between the continuum predictions and the numerical estimates. The situation changed when the 



6 But evidently, we rely on what we already learnt in previous subsections. 

7 It would be wonderful to find another way to impose the boundary, which preserves the BRST. Until now, we are 
unfortunately unaware of other consistent ways to do so. 



9 



The GZ framework 



David Dudal 



papers [39, Q appeared, where significantly larger lattices were invoked. Improved Schwinger- 
Dyson studies were able to reproduce these lattice results [|3(], [H], p2Q . But the Gribov-Zwanziger 
framework seemed to miss something essential. In [15, 41], we paid attention to other potentially 
important nonperturbative effects in the Gribov-Zwanziger theory, which could also influence the 
propagators. 

In particular, we focused on the dynamics of the additional fields, (Jp a ^ , <p„ c , ~c5 a £ ,(0^), which 
had to be introduced to consistently treat the theory. In a sense, these fields describe the influence 
of the Gribov horizon. Evidently, they will develop their own quantum dynamics, which can in 
return influence the conventional gluon and ghost sector, as these are nontrivially coupled to the 
"auxiliary" sector. In [ 15], we presented evidence for the existence of the dimension 2 condensate, 
ffiu < P'a ~ G^u 01 ??)- This operator is a very natural object to study, as the Gribov-Zwanziger theory 
already has an inherent mass scale in it, being the Gribov parameter y , and the horizon condition 
exactly corresponds to a d = 2 condensate, ( 1.19 ). Moreover, the operator is also compatible with 
the renormalizability [ ]4li |TJ]. To get a taste of the importance of such a condensate, we studied it 
in [15] using variational perturbation theory. Doing so, we found a gluon propagator of the type 

P 



D(p z 



p 4 +M 2 p 2 + X 4 '' 



while for the ghost, we obtained 



G(p 2 )~ forp 2 ~0. 



The mass parameter M 2 corresponds to the condensate (<p^ c <p^ c 
consistent way through the variational perturbative approach [H 



(6.1) 
(6.2) 

>, and it is fixed in a self- 



Both predictions of the so-called 
Refined Gribov-Zwanziger framework [ |15| ] are now in compliance with the lattice data. 



We can reexpress the gluon propagator (63) as 

1 



Dip 1 ) 



p 2 +m 2 (p 2 



with m (p 



I 4 



+ M 2 ' 



(6.3) 



This kind of "effective gluon mass" m (p ), which vanishes in the ultraviolet, has also been found 
in Schwinger-Dyson studies [ [30| , 12], and finds phenomenological application, see e.g. [43]. 

We are currently studying the issue of d = 2 condensates in the Gribov-Zwanziger theory 
more thoroughly, by using an effective potential approach for local composite operators [24]. We 
are not only including the aforementioned one, but also other ones like {(p^ c (p^ c ), see also [[25]]. 
If a nonvanishing vacuum expectation values for such condensates is dynamically favoured as the 
vacuum energy is lowered, then this would be an excellent illustration of the fact that the dynamics 
of the Gribov-Zwanziger theory is very rich, and needs to be taken into proper account to obtain 
reliable estimates for e.g. the propagators. Before finishing, it is perhaps important to point out that 
none of these condensates are free to choose mass parameters [31 1, each one of them is eventually 
expressed in terms of Aqcd by dimensional transmutation 8 . 

An issue we did not touch in this proceeding is the question what the physical operators might 
be in the Gribov-Zwanziger theory, since we do not longer have the BRST symmetry. We refer to 
p4 , 45, [l(|] for a discussion thereof. 

8 Notice that the horizon condition also fixes the Gribov parameter y 2 in terms of Aqcd |10[]- 
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